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Two families of models of rotating relativistic disks based on Taub- 
NUT and Kerr metrics are constructed using the well-known "dis- 
place, cut and reflect" method. We find that for disks built from 
a generic stationary axially symmetric metric the "sound velocity", 
(pressure /density) 1 ^ 2 , is equal to the geometric mean of the prograde 
and retrograde geodesic circular velocities of test particles moving on 
the disk. We also found that for generic disks we can have zones with 
heat flow. For the two families of models studied the boundaries that 
separate the zones with and without heat flow are not stable against 
radial perturbations (ring formation). 
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1 Introduction 



Axially symmetric solutions of Einstein field equations corresponding to disk 
like configurations of matter are of great astrophysical interest and have 
been extensively studied. These solutions can be static or stationary and 
with or without radial pressure. Solutions for static disks without radial 
pressure were first studied by Bonnor and Sackfield |IJ, and Morgan and 
Morgan ||, and with radial pressure by Morgan and Morgan ||, and, in 
connection with gravitational collapse, by Chamorro, Gregory and Stewart 
0]. Disks with radial tension has been recently considered in ||. Also models 
of disks with, electric fields 0, magnetic fields |7|], and both magnetic and 
electric fields have been introduced recently ||. Using a solitonic technique 
Neugebauer and Meinel found solutions representig rigidly rotating disks of 
dust H . Several classes of exact solutions of the Einstein field equations that 
represent disks were obtained by different authors [|10|1 - [|18[ . 



The stability of static models with no radial pressure can be explained 
by either assuming the existence of hoop stresses or that the particles on the 
disk plane move under the action of their own gravitational field in such a 
way that as many particles move clockwise as counterclockwise. This last 
interpretation is frequently made since it can be invoked to mimic true ro- 
tational effects. Even though, this interpretation can be seen as a device, 
there are observational evidence of disks made of streams of rotating and 



counter-rotating matter 19 



For disks without tension, i.e., usual fluid disks, the rotation is a necessary 
ingredient to have stability. Despite its relevance the literature on exact 
relavistic disks solutions to the Einstein field equations is scarce, we only 
found a discussion of a rotating disk as a source of Kerr metric PD . 



The aim of this work is to study models of rotating relativistic disks in 
some detail. The models are based in Taub-NUT and Kerr metrics, that are 
some of the simplest axially symmetric stationary solutions of the vacuum 
Einstein equations [J21| . To construct the disk solutions we use the well-known 
"displace, cut and reflect" method. 

We find that for disks built from a generic stationary axially symmetric 
metric the "sound velocity", (pressure/density) 1 ^ 2 , is equal to the geometric 
mean of the prograde and retrograde geodesic circular velocities of a test 
particle moving on the disk. We also found that the generic disk may have 
zones with heat flow. For the two families of models studied the boundaries 
that separates the zones with and without heat flow are not stable against 
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radial perturbations (ring formation). 

In Sec. II we study the generation of a disk model using the "displace, cut 
and reflect" method from a vacuum solution of the Einstein field equations 
for a stationary axially symmetric metric, in particular, we study: the general 
form of the associated energy-momentum tensor, disk velocity, and angular 
momentum. In the next section, Sec. Ill, we discuss the relation between 
the geodesic circular velocities on the disk and the pressure and density. In 
Sec. IV a class of models of rotating disks based on the Taub-NUT metric is 
presented. In particular, we study the apparition of heat flow and stability 
against radial perturbations. In Sec. V, another class of models of rotating 
disks is presented, this time based in the Kerr metric. Also the stability 
against radial perturbation is considered. The material presented here is 



complementary to the work of reference [20]. In the last section, Sec VI , we 



summarize our main results and made some comments. 



2 Rotating Relativistic Disks 

In this section we present a summary of the main quantities associated to 
the disk, we follow closely our reference ||. A sufficiently general metric for 
our purposes is Weyl-Lewis-Papapetrou line element, 

ds 2 = - e 2 "(dt + Ady) 2 + e" 2 W + e 2A (dr 2 + dz% (1) 

where $, A and A are functions of r and z only. Assuming the existence 
of the second derivatives of the functions A and A, the Einstein vacuum 
equations for this metric are equivalent to, 

^rr + ^ + ^ + ^i^ + Al) = 0, (2a) 
A 

A, rr + A ;ZZ + 4($, r ^ r + $, z A tZ ) = 0, (2b) 



„4$ 



A 



rK-^)-^{Al-Al), 



,4$ 



A i2 = 2r® >r $ }Z - — A, r A jZ . 



(2c) 
(2d) 
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Now if the first derivatives of the metric tensor are not continuous on the 
plane z = with discontinuity functions, 

bab = 9ab,z\ z=Q+ 9ab,z\ z=Q -i 

the Einstein equations yield an energy-momentum tensor (EMT) T ab = 
Qab fi(z), where 5(z) is the usual Dirac function with support on the disk 
and 

Qt = \{b az 5 z b - b zz 5 a b + g az ¥ b - g zz b a b + ¥ c (g zz 5 a b - g az 5 z b )}. 

is the distributional energy-momentum tensor. The "true" surface energy- 
momentum tensor of the disk can be written as S ab = e A ~* Q ab . 
For the metric ([I]) we obtain 

5 00 = 2e 3 *- A [2$, - AJ , (3a) 

5 01 = e 3 ^ A [4A<^ tZ -2AA, z + A, z ] (3b) 

S n = 2e 3 *- A [(rV 4 *-^ 2 )A, z + 2^ 2 $, + A4j (3c) 

where all the quantities are evaluated at z = + . 

The eigenvalue problem for the energy- momentum tensor ([3a|) - (|3c]), 

sz e = a e, (4) 

has the solutions 

A± = Ut±Vd), (5) 
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and A r = \ z = 0, where 



T = S° Q + Si D = (Si - S° ) 2 + 4 S? S 1 ^ (6) 
We can write 

g ab = -V a V b + W a W b + X a X b + Y a Y b , (7) 
and the canonical form of the EMT, 

Sab = aV a V b + PW a W b + K(V a W b + W a V b ), (8) 
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with a orthonormal basis 

V a = N o (l,n,0,0), (9a) 



(A- 


- Sg)/S? 


, D>0, 


(Sl~ 


- ^o )/2^ 


, D<0, 


(A+ 


-si)/ si 


D > , 




o 


D < , 



PF 1 = iV\(A, 1,0,0), (9b) 

X a = e *-A( 0; 0,1,0), (9c) 

Y a = e*~ A (0, 0,0,1), (9d) 

where Nq and JVi are normalization factors, and 

f (A_-5g)/S? , £ 

n = ^ (10) 



(11) 



The energy density, the azimuthal pressure, and the heat flow function 
are, respectively, 

-A_ , D > , 

(12) 

-T/2 , D < , 

A+ , D > , 
P = < (13) 
T/2 , £) < , 

, D > , 
= < (14) 
v^/2 , TJ < . 

The orthonormal basis {V a , W a , X a , Y a } is comoving with the disk. Thus 
the time-like vector V a will define the velocity vector of the disk: 

v a = (v°,v\o,o) = v°{i,n,o,o), (is) 
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where 

y° = . : en) 



e 



and 



V = , g "" 2 + 2gMii (17) 
V -^oo 

is the tangential velocity of the disk. 

The specific angular momentum of a particle of the disk, with mass /i, is 
given by 

u - P<p 



Thus, 



9 va V a . (18) 
A 4 



h = - 9uQ + 901 =. (19) 



A condition of stability under radial perturbations is 

^ = 2/^>0, (20) 
dr dr 

For h > 0, we have stability when the specific angular momentum is an 
increasing function of r. This criteria is an extension of Rayleigh criteria of 
stability of a fluid at rest in a gravitational field, see for instance p2| . 



3 The counter-rotating Model 

In this section we analyze the model of disks made of particles moving in 
opposite directions. When the discriminant D is positive there is not heat 
flow. Furthermore, when P > an observer comoving with the disks can 
consider the energy-momentum tensor as representing two streams of colli- 
sionless particles that circulate in opposite directions. 

Let be u a = (u°, u 1 , 0, 0) = u°(l, u, 0, 0) the velocity vector of the stream. 
The angular velocity uj can be obtained from the geodesic equation for a test 
particle, we get 

gu^ 2 + 2g 01iT uj + g 00jr = 0, (21) 
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with solutions, 



—g01,r ± \l 9oi,r ~ 900,rgil,r 
9ll,r 



(22) 



Therefore, in general, the two streams circulate with different velocities. 

We can compute the tangential velocity of the streams by projecting the 
velocity vector u a onto the comoving tetrad, ea 6 = {V b ,W b ,X b ,Y b }, 



u 



e\u b = V M e db u b . 



We get, 







W + WiUJ± 


vP 




V + V x w± 



By performing the product of U+ and LL we obtain 

W 2 + W W 1 (cu + + w_) + Wx 



U+U_ 



V 2 + V V 1 (uj + + uj_) + V^u+uj- 
From (|2ip we have the relations, 

2goi,r goo.r 



L0 + + LU_ 



9u, 



9n, 



Thus, 



g 11:r W 2 - 2g 01 -W W 1 + g^W^ 



9ii,rV 2 - 2g 0hr V V 1 + g m , r Vi 2 



By using fl7|) and (§), we can write 



A + aB 



A-PB 



where 



A — gn,rSoo — 2g i r Soi + goo, r Sn, 
B = goo9n,r — 2goi9oi,r + 900,r9ll- 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 

(29a) 
(29b) 
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Using the Einstein equations ( Pa]) - ( |2"dD and the expressions 
the energy-momentum tensor (with K = 0) we can show that 

A = -2r(S° + Sl). 

Also, using the metric (|l]) we get, 

B = -2r. 

Hence, assuming that P > and a > 0, we obtain 



XJJJ_ 



Si + Sl + a 



5*0° + SI-P 



P 



a 



( |3"c| ) for 

(30) 

(31) 

(32) 
(33) 



In other words, for disks built from a generic stationary axially symmetric 
metrics the "sound velocity" (P/a) 1 / 2 is equal to the geometric mean of the 
prograde and retrograde geodesic circular velocities of a test particle moving 
on the disk. 



4 Taub-NUT Disks 

The simplest stationary axially symmetric solution of the Einstein equations 
is the Taub-NUT solution that can be written as ([!]) with 

x 2 -l 



A 



A 



iln 

2 

bn 
2 



x 2 + 2px + 1 



x 



x 2 — y 2 



2aqy, 



(34a) 

(34b) 
(34c) 



where p = m/a , q = I /a , with a 2 = m 2 + / 2 , so that p 2 + q 2 = 1. m and / 
are the mass and the NUT parameter, respectively, x and y are the prolate 
spheroidal coordinates, given by 



2ax 



y/r 2 + (z + b + 



a) 



r 2 + (z + b — a) 2 , 



(35a) 



2ay 



yV 2 + (z + b + a) 2 - \/r 2 + (z + b- a)' 



(35b) 
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with b > a > 0. Note that we have displaced the origin of the z- axis in b. 
From the above expressions we can compute the trace (T) and discriminant 
(D) of the energy-momentum tensor. We found T = aT and D = a 2 D, with 



4y[y 2 (2x 3 + 3px 2 — p) + x(px 3 — 3px — 2)] 
\(x 2 -y 2 ){x 2 + 2px + l)] 3 / 2 



x and y are given by 



(36a) 



£ = 16[y 2 (px 2 + 2x + p) 2 - q 2 x\x 2 - 1)(1 - y 2 )] (3gb) 
(x 2 — y 2 )(x 2 + 2px + l) 3 



2x = 2 + (k + l) 2 + ^ 2 + (« - I) 2 , (37a) 



2y = + ( K + 1)2 - y/jx + ( K - 1)2, (37b) 

where r = r/a and k = bja. 

Even though we have an exact solution and all the relevant quantities of 
the disk can be explicitly computed the resulting expressions are cumber- 
some and not very illuminating, the adoption of the graphic method is more 
adequate in this case. The first quantity to be considered is the discriminant 
D that will give us the canonical form of the EMT. 

In Fig. [I] we show D for Taub-NUT disks with different values of p and 
k. First we plot D for disks with k = 1.7 and p = 0.9, 0.7, 0.5, 0.3, and 0.1. 
Then we plot D for disks with p = 0.8 and k — 2, 2.6, 3.2, 3.8, and 4.4. We 
also computed D for many other values of p and k. We found that in all the 
cases when p ^ 1, D is not positive definite and it has only one root f > 0. 
Therefore the Taub-NUT disks always have heat flow beginning at f = f . 

In Fig H we depict the density a = aa and the pressure P = aP for 
Taub-NUT disks with different values of p and k. First we plot a and P 
(scaled by a factor 10) for Taub-NUT disks with k = 1.4 and p = 1, 0.8, and 
0.6 as functions of f = r/m. Then we plot a and P (scaled by a factor 10) 
for Taub-NUT disks with k = 1.7 and p — 1, 0.8, and 0.6 as functions of 
r = r/m. We also computed a and P for many other values of p and k, in all 
the cases, we obtain a similar behaviour. The density a is always positive, 
a > 0; whereas the pressure becomes negative (tension) for a value of f < fo. 
We can also see that a and P presents a non smooth behaviour at f = f . 
For f > f we have heat flow. 
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The heat flow function K = aK is represented in Fig. []. We plot K for 
k = 1.4 (upper two curves) and k = 1.7 (lower two curves), with p = 0.8 and 
0.6 as functions of f. Also in Fig. |3], in order to see the change of behaviour 
of a and P at f = fo and the relation between a and K, we plot a, P and K 
for n = 1.1 and p = 0.8 in the interval 2.5 < f < 3. We can see that K > a 
for f > fi > f . Thus, there is not causal propagation of heat for f > fi. 

In Fig. f| we show the tangential disk velocity V and the disk angular 
momentum h for Taub-NUT disks with p = 0.8 and k = 1.1, 1.4, 1.7 and 2, 
as functions of f. We see that there is a strong change in the slope of V and 
h at f = fo- That means that there is a strong instability at this value of f. 
We also show the stream angular momenta h + and 7t_ in order to compare 
with the counter-rotating model. We also found instability, but in a different 
place. We computed also V and h for a wide range of the parameters p and 
k, we found always the same behaviour. 

In order to compare the counter-rotating model of the disk with the true 
disk rotation, we plot in Fig. |5] the tangential velocities of the counter- 
rotating streams, U + (top) and £/_, and the product U + U- (bottom full line) 
and the P/a (full dots) for k = 1.1 and p = 0.8, as functions of f. We can see 
the exact matching of the two curves, U + U- and P/a, in perfect agreement 
with (|33|) . Also in Fig. [|, in order to show the behavior of the counter- 
rotating model for the different disks models, we plot the P/a relation for 
k = 1.4 (upper curves) and k = 1.7 (lower curves), with p = 1, 0.8, and 0.6 
as functions of f. 



5 Kerr Disks 



In this section we use the "displace, cut, and reflect" method to built disk 
solutions using the Kerr solution, that can be written as ([!]) with 



A 



A 



iln 

2 



-In 

2 

2aq 
p 



p l x + q y — 1 
(px + l) 2 + q 2 y 2 

p 2 x 2 + q 2 y 2 — 1 
p 2 {x 2 — y 2 ) 

(1 — y 2 ){px + 1 



p 2 x 2 



q 2 y 2 



(38a) 
(38b) 
(38c) 
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where p = a/m , q = a/m , with a 2 = m 2 — a 2 , so that p 2 + q 2 = 1. m and 
a are the mass and the Kerr parameter, respectively, x e y are, again, the 
prolate spheroidal coordinates, given by (|35aj ) and ( |35b| ). 

With the above expressions we can compute the trace and the discrimi- 
nant of the energy-momentum tensor, that now can be written as T = mT 
and D = m 2 D, with 

f = 4xy(x 2 -l)(l-y 2 ) 

p{x 2 — y 2 ) 2 (p 2 x 2 + q 2 y 2 — 1) 



2q 2 x(px + 1)(1 — y 2 ) + p(x 2 — l)[(px + l) 2 — q 2 y 2 ] 
p{x 2 — y 2 ){p 2 x 2 + q 2 y 2 — l)[(px + l) 2 + q 2 y 2 ] 



jj = 4 - 2 / 2q 2 x(px + !)(!- y 2 ) + p(x 2 - 1) [(px + l) 2 - q 2 y 2 } 
\ p(x 2 — y 2 )(p 2 x 2 + q 2 y 2 — + l) 2 + q 2 y 2 ] 

16q 2 y 2 (x 2 - 1)(1 - y 2 )[3p 2 x 2 + Apx - q 2 y 2 + l] 2 



(39b) 



p 2 (x 2 — y 2 ) 2 (p 2 x 2 + q 2 y 2 — l) 3 [(px + l) 2 + q 2 y 2 } ' 
where x and y are given by 

2px = \/f 2 + (k + p) 2 + ^/f 2 + (k — p) 2 , (40a) 



2py = a/?- 2 + (k + p) 2 - >/r 2 + (k -p) 2 , (40b) 

where r = r/m and n = b/m. 

In Fig. ^ we show D for Kerr disks with different values of p and k. We 
first plot D for disks with k = 1.5 and p = 1, 0.9, 0.8, 0.6, and 0.3. Then 
we show D for disks with p = 0.8 and k = 1.7, 2, 2.3, 2.9, and 3.8. We also 
computed D for different values of p and k, we find that for some values of p 
and k, D > everywhere on the disk, whereas for the complementary case 
D have two roots, < r\ < f 2 , so that D < when < fx < f 2 . Accordingly 
we have two kinds of Kerr disks, ones with heat flow in an annular region 
(0 < fx < f 2 ) and others without heat flow. 

We first analyze some cases when there is no heat flow. In Fig. [F] we show 
the density a = ma and the pressure P = mP for these Kerr disks. First we 
plot a and P (scaled by a factor 100) for a disk with p = 0.9 and k = 2.4, 
2.5, and 3 as functions of f = r/m. Then we plot a and P (scaled by a factor 
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10) for disks with k = 2.5 and p = 1, 0.9, and 0.8. We also computed a and 
P for a variety of p and k, we found a similar behaviour. The density a is 
positive definite everywhere on the disk, falling to zero at infinity. There are 
values of p and k such that the pressure P is positive definite, also falling to 
zero at infinity, whereas in other cases the pressure P is negative (tension) 
at the central region of the disks becoming positive for a value f > 0; then 
falling to zero at infinity. 

The rotational motion of these disks models is shown in Fig. §. We first 
plot V for p = 0.9 with k = 2.4, 2.5, and 3 and for p = 0.8 with k = 2.5 (top 
curve) and then we plot the disk angular momentum h for the same values 
of p and k. We can see that the disk velocity increases initially and then falls 
to zero at infinity and always is less than the light velocity. Also we see that 
the disk angular momentum is an increasing monotonic function of r. Thus 
these disks models are stable. 

Finally, we analyze a Kerr disk with heat flow, with p = 0.9 and n = 1.5. 
In Fig. ^ we show a, P and K for this values of p and k. We see that the 
density a is positive definite everywhere on the disk, falling to zero at infinity 
and that presents a change of slope at the values of r when D — 0, f\ and 
T2- The pressure P is negative at the central region of the disk, becoming 
positive for a value of f > T2- Also, P presents a change of slope at f\ and fi- 
The heat flow function K increases rapidly and becomes larger than a for a 
value of f > fi, K falls to zero, becoming less than a for a value of r < f 2 . 
So, there is not causal propagation of heat on this region of the disk. 

We also plot in Fig. |9] the disk angular momentum, h. We see that 
the angular momentum presents two instabilities, at the values of f when 
D — 0, fi and f2, and so the region of the disk where there is heat flow is 
highly instable. In order to see that these instabilities are present also in the 
counter-rotating model, we plot also the stream angular momenta, h + and 
/i_ (scaled by a factor 20). We analyzed several other cases of Kerr disks 
with heat flow, we found the same behaviour. 

6 Discussion 

We presented two families of rotating disks based on the Taub-NUT and 
Kerr solutions. For the first family we find that the presence of heat flow is a 
generic property of this models, also generic is the radial instability produced 
by this flow. For the second family we find two sub classes one with and the 
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other without heat flow. The heat flow is concentrated in an annular region 
that is highly instable. 

The relation found between the geodesic circular velocities and the "speed 
of sound" found is a generalization for the stationary case of the correspond- 
ing relation found by Morgan and Morgan || that is the base of the counter- 
rotating model. 

The inclusion of current in the rotating models to have a rotating "hot" 
disk as the ones studied in |7[ is under consideration. The rotation add mass 
to the otherwise massless charges. These models are quite nontrivial. Also, 
the inclusion of radial pressure or tension is being considered. 

Finally, we want to mention that all the computation of this work was 
performed using the algebraic programming system Reduce 153 . 
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Figure 1: We plot D for Taub-NUT disks with « = 1.7 and p = 0.9, 0.7, 0.5, 
0.3, and 0.1 as functions of f = r/a. Next we plot D for Taub-NUT disks 
with p = 0.8 and k — 2, 2.6, 3.2, 3.8, and 4.4 as functions of f = r/a. Note 
that -D is not positive definite and has only one root f > 




Figure 2: First we plot a and P (scaled by a factor 10) for Taub-NUT disks 
with k — 1.4 and p—1, 0.8, and 0.6 as functions oi r = r /m. Next we plot cr 
and P (scaled by a factor 10) for a Taub-NUT disk with k — 1.7 and p — 1, 
0.8, and 0.6 as functions of f = r/m. P presents a non smooth behaviour at 
f — f , wherein the heat flow begins. 
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Figure 3: We show K for k = 1.4 (upper curves) and k — 1.7 (lower curves), 
with p = 0.8 and 0.6 as functions of f. Next we plot a, P and for k — 1.1 
and p = 0.8 in the interval 2.5 < f < 3. 




Figure 4: We plot the tangential disk velocity, V, with p = 0.8 and k — 1.1, 
1.4, 1.7, and 2 as functions of f. Next we plot, for the same values of p and 
k , the disk angular momentum, —h (sharp curves), and the stream angular 
momenta, h+ (top curve from the left) and — h- (scaled by a factor of 15), 
for p = 0.8 and k = 1.1 as functions of f. 
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U+, U-, U+U-, P/a P/a 




r = r/a f — rj 



Figure 5: We show U + , LL, U + U- and P/a (full dots) for k = 1.1 andp = 0.8 
as functions of f. Next we plot the P/a relation for k = 1.4 (upper curves) 
and k — 1.7 (lower curves), with p—1, 0.8, and 0.6 as functions of f. 




Figure 6: The discriminant D associated with the Kerr disks with k = 1.5 
is plotted for p — 1, 0.9, 0.8, 0.6, and 0.3. Next we show D for disks with 
p = 0.8 and k = 1.7, 2, 2.3, 2.9 and 3.8. Note that D is not positive defined. 
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Figure 7: First we plot a and P (scaled by a factor of 100) for Kerr disks 
with p = 0.9 and k = 2.4, 2.5, and 3 as functions of f — r / m. Next we plot 
a and P (scaled by a factor of 10) for disks with k = 2.5 and p = I, 0.9 and 
0.8. 




Figure 8: We plot V for p = 0.8 with k = 2.5 (top curve) and p = 0.9 with 
k = 2.4, 2. 5, and 3 (following three curves). Next we plot the disk angular 
momentum, h. for the same values of p and k. 
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Figure 9: We plot the density , a, the pressure, P, and the heat function 
function, K, for a Kerr disk withp = 0.9 and k = 1.5 as functions of f = r/m. 
Next, for the same values of p and k, we plot the disk angular momentum, 
h, and the stream angular momenta, h + and h- (scaled by a factor 20). 
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